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This paper is devoted to a systematic exposure of the Kondo physics in quantum dots for which 
the low energy spin excitations consist of a few different spin multiplets \SiMj). Under certain 
conditions (to be explained below) some of the lowest energy levels Es i are nearly degenerate. 
The dot in its ground state cannot then be regarded as a simple quantum top in the sense that 
beside its spin operator other dot (vector) operators R n are needed (in order to fully determine its 
quantum states), which have non-zero matrix elements between states of different spin multiplets 
(SiMi\R, n \SjMj) 7^ 0. These "Runge-Lenz" operators do not appear in the isolated dot-Hamiltonian 
(so in some sense they are "hidden"). Yet, they are exposed when tunneling between dot and leads 
is switched on. The effective spin Hamiltonian which couples the metallic electron spin s with the 
operators of the dot then contains new exchange terms, J n s ■ R n beside the ubiquitous ones JiS ■ Si. 
The operators Si and R n generate a dynamical group (usually SO(n)). Remarkably, the value of n 
can be controlled by gate voltages, indicating that abstract concepts such as dynamical symmetry 
groups are experimentally realizable. Moreover, when an external magnetic field is applied then, 
under favorable circumstances, the exchange interaction involves solely the Runge-Lenz operators 
R n and the corresponding dynamical symmetry group is SU(n). For example, the celebrated group 
SU(3) is realized in triple quantum dot with four electrons. 

PACS numbers: 72.10.-d, 72.15.-v, 73.63.-b 



I. INTRODUCTION 

Recently, studies of the physical properties of artifi- 
cially fabricated nano-objects turn out to be a rapidly 
developing branch of fundamental and applied physics. 
Progress in these fields is stimulated both by the achieve- 
ments of nanotechnology and by the ambitious projects of 
information processing, data storage, molecular electron- 
ics and spintronics. The corresponding technological evo- 
lution enabled the fabrication of various low-dimensional 
systems from semiconductor heterostructures to quan- 
tum wires and constrictions, quantum dots (QD), molec- 
ular bridges and artificial structures with large molecules 
built in electric circuitsi. This impressive experimental 
progress led to the development of nanophysics, a new as- 
pect and research direction in quantum physics^. Artifi- 
cial nano-objects possess the familiar features of quantum 
mechanical systems, but sometimes one may create in 
artificially fabricated systems such conditions, which are 
hardly observable "in natura". For example, ID — > 2D 
crossover may be realized in quantum networks^ and 
constrictions^. The Kondo effect may be observed in non- 
equilibrium conditions^, at high magnetic fields^, and 
at finite frequencies^. Moreover, a quantum dot in the 
Kondo regime can be integrated into a circuit exhibiting 
the Aharonov-Bohm effect^. 

In this paper we focus on an intriguing challenge in 
this context related to the specific symmetry of the nano- 
objects under study. More precisely, one is interested in 
answering questions pertaining to the nature of the un- 
derlying symmetry of the dot Hamiltonian and the alge- 
bra of operators appearing in the exchange Hamiltonian. 
The investigation of this topic is intimately related with 
the geometric structure and electron occupation of the 



quantum dot in its ground state. We refer to a quan- 
tum dot composed of a single well and containing an odd 
number N of electrons as a simple quantum dot (SQD). 
The dot Hamiltonian of a SQD ( in the absence of an 
external magnetic field) is composed of two degenerate 
levels and has an SU(2) symmetry. In that sense the 
symmetry is referred to as geometrical. The exchange 
Hamiltonian is expressed in terms of the generators of 
the group SU(2). On the other hand, quantum dots con- 
taining a single well with even N or quantum dots con- 
taining several wells are referred to as composite quantum 
dots (CQD). The low energy states of an isolated CQD 
Hamiltonian are spin multiplets. In the generic case, the 
only degeneracy is that of magnetic quantum numbers. 
Yet, as we argue below, dot-lead tunneling results in level 
renormalization and the emergence of an additional de- 
generacy, both generic and accidental. To be more pre- 
cise, we note that 1) The exchange part of the Hamilto- 
nian includes the generators of a non-compact Lie group 
(usually SO(n) or SU(n)) and 2) The renormalized low- 
energy spin-excitation levels of the CQD Hamiltonian are 
almost degenerate (within a Kondo energy scale). These 
two aspects are gathered under the term dynamical sym- 
metry. A more quantitative exposition will be presented 
below. 

Experimentally, resonance Kondo tunneling was ob- 
served in QD with odd electron occupation number under 
strong Coulomb blockade^, and in individual atoms and 
molecules deposited on metallic surfaces and on the edges 
of metallic wires in break-junction geometry^. Accord- 
ing to the theory of Kondo effect in QD 11 , spin degrees 
of QD are involved in Kondo resonance. In our notation 
these are SQD and the physics of Kondo tunneling in 
this case is similar to that of Kondo scattering in mag- 
netically doped metals, at least in the regime of linear 
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response. 

The Kondo physics seems to be richer in systems in- 
volving tunneling through CQD. Our main purpose is to 
demonstrate that CQD possesses dynamical symmetries 
whose realization in Kondo tunneling is experimentally 
tangible. Such experimental tuning of dynamical sym- 
metries is not possible in conventional Kondo scattering. 
In many cases even the very existence of Kondo tunnel- 
ing crucially depends on the dynamical spin symmetry 
of CQD. Several models dealing with Kondo tunneling 
in CQDs possessing dynamical symmetry were consid- 
ered in our previous publications. The case of 5*0(4) 
symmetry in double quantum dot (DQD) was studied in 
Refs. ll2l[T^ A more complicated case of SO(n) sym- 
metry with variable n in a triple quantum dot (TQD) 
(composed of three potential valleys with intra-well in- 
teraction and inter-well tunneling) is introduced in Ref. 

M 

The main goal here is to develop the general approach 
to the problem of dynamical symmetries in Kondo tun- 
neling through CQDs and illustrate it by numerous ex- 
amples of TQD in various configurations, both in parallel 
and in series geometries. Within this framework, our ear- 
lier results^ fit into an elegant pattern of classification of 
dynamical symmetry groups, which is expanded here in 
a somewhat more complete and rigorous formalism. The 
main lesson to be learned is that Kondo physics in CQD 
suggests a novel and in some sense rather appealing as- 
pect of low-dimensional physics of interacting electrons. 
It substantiates, in a systematic way, that dynamical 
symmetry groups play an important role in mesoscopic 
physics. In particular, we encounter here some "famous" 
groups which appear in other branches of physics. Thus, 
the celebrated group 5t/(3) enters also here when a TQD 
is subject to an external magnetic field. And the group 
50(5) which plays a role in the theory of superconduc- 
tivity is found here when a certain tuning of the gate 
voltages in TQD is exercised. 

The basic concepts are introduced in section [H] First, 
in sub-section 111 Al the necessary mathematical ingredi- 
ents are introduced, although we try to avoid much rigor. 
Then, in sub-section III Bl we explain how these abstract 
concepts can be realized in CQD. In section ITTTl the spe- 
cial case of TQD in the parallel geometry is discussed 
at some length. In subsection IIII Al we derive scaling 
equations for TQD with even occupation and calculate 
Kondo temperatures for various dynamical symmetries. 
In section HV1 we discuss the physics of TQD in a series 
geometry and point out similarities and differences be- 
tween Kondo physics in the two geometries (sub-section 
IIV A|) . In sub-section IIVBI we concentrate on the case 
of even occupation. The dynamical-symmetry phase dia- 
gram is displayed and the experimental consequences are 
drawn. The case of odd occupation is exposed in sub- 
section IIV CI Finally, in section a novel phenomena 
is discussed, namely, a Kondo effect without a localized 
spin. The anisotropic exchange interaction occurs be- 
tween the metal electron spin and the dot Runge-Lenz 



operator alone. In the conclusions we underscore the 
main results obtained here. 

The derivation of the pertinent effective spin Hamiltoni- 
ans and the establishment of group properties (in particu- 
lar identification of the group generators and checking the 
corresponding commutation relations) sometimes require 
lengthy mathematical expressions. These are collected in 
the appendices. 



II. DYNAMICAL SYMMETRY OF COMPLEX 
QUANTUM DOTS 

A. Generalities 

In this section we present in some details the concept 
of dynamical symmetry, and more particularly, its emer- 
gence in CQD. The term Dynamical Symmetry implies 
the symmetry of eigenvectors of a quantum system form- 
ing an irreducible representations of a certain Lie group. 
The main ideas and the relevant mathematical tools can 
be found, e.g., in Refs. 0. Here they are formulated in a 
form convenient for our specific purposes without much 
mathematical rigor. We have in mind a quantum system 
with Hamiltonian H whose eigenstates |A) = \Mfi) form 
( for a given M) a basis to an irreducible representation 
of some Lie group Q. The energies Em do not depend 
on the "magnetic" quantum number \i. For definiteness 
one may think of M as an angular momentum and of /i 
as its projection, so that Q is just SU(2). Now let us 
look for operators which induce transitions between dif- 
ferent eigenstates. An economic way for identifying them 
is through the Hubbard operators^ 

X AA ' = |A)(A'|. (1) 

It is natural to divide this set of operators into two sub- 
sets . The first one contains the operators \M fi) (fjf M\ 
while the second one includes operators \Mfj,)(fj/M'\ for 
which |M/i) and \M'fi'} belong to a different represen- 
tation space of Q. A central question at this stage is 
whether these operators (or rather, certain linear combi- 
nations of them) form a close algebra. In some particu- 
lar cases it is possible to form linear combinations within 
each set and obtain two new sets of operators {5} and 
{R} with the following properties: 1) For a given M the 
operators {5} generate the M irreducible representation 
of Q and commute with H. 2) For a given set Mi the 
operators {5} and {R} form an algebra (the dynamic al- 
gebra) and generate a non-compact Lie group A. The 
reason for the adjective dynamic is that, originally, the 
operators {R} do not appear in the bare Hamiltonian H 
and emerge only when additional interaction (e.g., dot- 
lead tunneling) is present. In the special case Q = SU{2) 
the operators in {5} are the vector S of spin operators de- 
termining the corresponding irreducible representations, 
while the operators in the set {R} can be grouped into 
a sequence R„ of vector operators describing transitions 
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between states belonging to different representations of 
the SU(2) group. 

Strictly speaking, the group A is not a symmetry group 
of the Hamiltonian H since the operators {R} do not 
commute with H. Indeed, let us express H in terms of 
diagonal Hubbard operators, 

H= ]T £a|A}(A|=5> m X AA , (2) 

A=iW> A 

so that 

[X AA ' ,H} = - (E M - E M - )X aa ' . (3) 

As we have mentioned above, the symmetry group Q 
of the Hamiltonian H, is generated by the operators 
X A=M ^' A = m m Remarkably, however, the dynamics of 
CQD in contact with metallic leads and/or an external 
magnetic field leads to renormalization of the energies 
{Em} in such a way that a few levels at the bottom of the 
spectrum become degenerate, Em x = Em 2 = • • • Em„ ■ 
Hence, in this low energy subspace, the group A gener- 
ated by the operators {S} and {R} is a symmetry group 
of H referred to as the dynamical symmetry group. The 
symbol R is due to the analogy with the Runge-Lenz op- 
erator, the hallmark of dynamical symmetry of the Ke- 
pler and Coulomb problems 17 . Below we will use the 
term dynamical symmetry also in cases where the levels 
are not strictly degenerate, their difference is bounded by 
a certain energy scale, which is the Kondo energy in our 
special case. In that sense, the symmetry is of course not 
exact, but rather approximate. 

Using the notions of dynamical symmetry, numerous 
familiar quantum objects, such as hydrogen atom, quan- 
tum oscillator in d-dimensions, quantum rotator, may be 
described in a compact and elegant way. We are inter- 
ested in a special application of this theory, when the 
symmetry of the quantum system is approximate and its 
violation may be treated as a perturbation. This aspect 
of dynamical symmetry was first introduced in particle 
physics 18 , where the classification of hadron eigenstates 
is given in terms of non-compact Lie groups. In our case, 
the rotationally invariant object is an isolated quantum 
dot, whose spin symmetry is violated by electron tunnel- 
ing to and from the leads under the condition of strong 
Coulomb blockade. 



B. Realization in CQD 

The special cases Q = SU(2) and A = SO(n) or SU(n) 
is realizable in CQD. Let us first recall the manner in 
which the spin vectors appear in the effective low en- 
ergy Hamiltonian of the QD in tunneling contact with 
metallic leads. When strong Coulomb blockade com- 
pletely suppresses charge fluctuations in QD, only spin 
degrees of freedom are involved in tunneling via the 
Kondo mechanismii. An isolated SQD in this regime 



is represented solely by its spin vector S. This is a man- 
ifestation of rotational symmetry which is of geometrical 
origin. The exchange interaction Js • S (s is the spin oper- 
ator of the metallic electrons) induces transitions between 
states belonging to the same spin (and breaks SU (2) in- 
variance). On the other hand, the low energy spectrum 
of spin excitations in CQD is not characterized solely 
by its spin operator since there are states close in en- 
ergy, which belong to different representation spaces of 
SU(2). Incidentally, these might have either the same 
spin S (like, e.g, in two different doublets) or a different 
spin (like, e.g., in the case of singlet-triplet transitions). 
The exchange interaction must then contain also other 
operators R„ (the R-operators mentioned in the previ- 
ous sections) inducing transitions between states belong- 
ing to different representations. The interesting physics 
occurs when the operators R„ "approximately" commute 
with the Hamiltonian Hdot of the isolated dot. In accor- 
dance with our previous discussion, the R-operators are 
expressible in terms of Hubbard operators and have only 
non-diagonal matrix elements in the basis of the eigen- 
states of Hdot. The spin algebra is then a subalgebra of 
a more general non-compact Lie algebra formed by the 
whole set of vector operators {S,R„}. This algebra is 
characterized by the commutation relations, 

[Si,Sj] = itijkSk, [Si, Rnj] = itijkRnki 

[R n i,R n j] = itj^jf.Sk, (4) 

with structure constants tijk, tfj k (here ijk are Cartesian 
indices). The R-operators are orthogonal to S, 

S • R„ = 0. (5) 

In the general case, CQDs possess also other symme- 
try elements (permutations, reflections, finite rotations). 
Then, additional scalar generators A p arise. These gener- 
ators also may be expressed via the bare Hubbard oper- 
ators, and their commutation relations with R-operators 
have the form 

if™ p R mj , (6) 

with structure constants g™ mp and /™ mp [n ^ to). The 
operators obeying the commutation relations I0J and © 
form an o„ algebra. The Casimir operator for this algebra 
is 

/C = S 2 + 5>2+£^. (7) 

n p 

Various representations of all these operators via basic 
Hubbard operators will be established in the following 
sections, where the properties of specific CQDs are stud- 
ied. 

Next, we show how the dynamical symmetry of CQD 
is revealed in the effective spin Hamiltonian describing 
Kondo tunneling. This Hamiltonian is derived from the 
generalized Anderson Hamiltonian 

Ha = Hdot + Hi ea d + H tU n ■ (8) 
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The three terms on the right hand side are the dot, lead 
and tunneling Hamiltonians, respectively. In the generic 
case, a planar CQD is a confined region of a semicon- 
ductor, with complicated multivallcy structure secluded 
between drain and source leads. The CQD contains sev- 
eral valleys numbered by index a. Some of these valleys 
are connected with each other by tunnel channels charac- 
terized by coupling constants W aa > , and some of them are 
connected with the leads by tunneling. The correspond- 
ing tunneling matrix elements are V ao (b = s,d stands 
for source and drain, respectively). The total number of 
electrons N in a neutral CQD as well as the partial occu- 
pation numbers N a for the separate wells are regulated 
by Coulomb blockade and gate voltages v ga applied to 
these wells, with N = ^2„N a . It is assumed that the 
capacitive energy for the whole CQD is strong enough 
to suppress charged states with N' = N ± 1, which may 
arise in a process of lead-dot tunneling. 

If the inter- well tunnel matrix elements W aa > are larger 
than the dot-lead ones V a b (or if all tunneling strengths 
are comparable), it is convenient first to diagonalize Hdot 
and then consider H tun as a perturbation. In this case 
Hdot may be represented as 



Thus, 



H dot = ^a|A)(A| + ]T E X \X)(X\. 



(9) 



AeN 



\£N±1 



Here all intradot interactions are taken into account. The 
kets | A) = \N,q) represent eigenstates of H dot in the 
charge sector N and quantum numbers q, whereas the 
kets | A) = \N± l,p) are eigenstates in the charge sectors 
N ± 1 with quantum numbers p. All other charge states 
are suppressed by Coulomb blockade. Usually, q and p 
refer to spin quantum numbers but sometimes other spec- 
ifications are required (see below). 

The lead Hamiltonian takes a form 



E 



(10) 



In the general case, the individual dots composing the 
CQD are spatially separated, so one should envisage the 
situation when each dot is coupled by its own channel to 
the lead electron states. So, the electrons in the leads 
are characterized by the index a, which specifies the lead 
(source or drain) and the tunneling channel, as well as 
by the wave vector k and spin projection a. 

The tunnel Hamiltonian involves electron transfer be- 
tween the leads and the CQD, and thus couples states 
| A) of the dot with occupation N and states |A) of the 
dot with occupation N ± 1. This is best encoded in terms 
of non-diagonal dot Hubbard operators, which intermix 
the states from different charge sectors 



ff * = E E (v a i X cl ka \A)(X\+H, 

+ E E (K» A M fe JA)<A|+i7.c.),(12) 



kaaa \£N-1,A£N 

where = V a {X\d aa \K). 

Before turning to calculation of CQD conductance, the 
relevant energy scales should be specified. First, we sup- 
pose that the bandwidth of the continuum states in the 
leads, D a , substantially exceeds the tunnel coupling con- 
stants, D a >• W aa ' j V a (actually, we consider leads made 
of the same material with D as — D a d — -Do)- Second, 
each well a in the CQD is characterized by the "activa- 
tion energy" defined as A a = EA(N a ) — E\(N a — 1), i.e., 
the energy necessary to extract one electron from the well 
containing N a electrons and move it to the Fermi level 
of the leads (the Fermi energy is used as the reference 
zero energy level from now on). Note that A a is tunable 
by applying the corresponding gate voltage v ga - We are 
mainly interested in situations where the condition 



A, ~ D () 



(13) 



is satisfied at least for one well labeled by the index c. 
Here Q c is a capacitive energy, which is predetermined 
by the radius of the well c. Eventually, this well with 
the largest charging energy is responsible for Kondo-likc 
effects in tunneling, provided the occupation number N c 
is odd. The third condition assumed in our model is a 
weak enough Coulomb blockade in all other wells except 
that with a = c, i.e., Q a <C Q c - Finally, we demand that 



b aa = ^-<l, 



(14) 



for those wells, which are coupled with metallic leads, 
and 



Pa = 



E„, 



< 1. 



(15) 



X 



AA 



|A)(A|, X 



AA 



|A)(A|. 



(11) 



Here E ac are the charge transfer energies for electron 
tunneling from the c-well to other wells in the CQD. 

The interdot coupling under Coulomb blockade in each 
well generates indirect exchange interactions between 
electrons occupying different wells. Diagonalizing the dot 
Hamiltonian for a given N = Y] N a , one easily finds 
that the low-lying spin spectrum in the charge sectors 
with even occupation N consists of singlet /triplet pairs 
(spin S = or 1, respectively). In charge sectors with 
odd N the manifold of spin states consists of doublets 
and quartets (spin S—l/2 and 3/2, respectively). 

The resonance Kondo tunneling is observed as a 
temperature dependent zero bias anomaly in tunnel 
conductance^. According to existing theoretical under- 
standing, the quasielastic cotunneling accompanied by 
the spin flip transitions in a quantum dot is responsible 



5 



for this anomaly. To describe the cotunneling through 
a neutral CQD with given N, one should integrate out 
transitions involving high-energy states from charge sec- 
tors with N' = N ± 1. In the weak coupling regime at 
T > Tk this procedure is done by means of perturbation 
theory which can be employed in a compact form within 
the renormalization group (RG) approach formulated in 
Refs. llHEl 

As a result of the RG iteration procedure, the energy 
levels E\ in the Hamiltonian @ are renormalized and 
indirect exchange interactions between the CQD and the 
leads arise. The RG procedure is equivalent to summa- 
tion of the perturbation series at T > Tk, where Tk 
is the Kondo energy characterizing the crossover from a 
perturbative weak coupling limit to a non-perturbative 
strong coupling regime. The leading logarithmic approx- 
imation of perturbation theory corresponds to a single- 
loop approximation of RG theory. Within this accuracy 
the tunnel constants W and V are not renormalized, as 
well as the charge transfer energy A c Reduction 
of the energy scale from the initial value Dq to a lower 
scale ~ T results in renormalization of the energy lev- 
els E\ — ► E\(Dq/T) and generates an indirect exchange 
interaction between the dot and the leads with an (anti- 
ferromagnetic) exchange constant J. 

The rotational symmetry of a simple quantum dot is 
broken by the spin-dependent interaction with the leads, 
which arises in second order in the tunneling amplitude 
V a . In complete analogy, the dynamical symmetry of a 
composite quantum dot is exposed (broken) as encoded 
in the effective exchange Hamiltonian. In a generic case, 
there are, in fact, several exchange constants arranged 
within an exchange matrix J which is non-diagonal both 
in dot and lead quantum numbers. The corresponding 
exchange Hamiltonian is responsible for spin-flip assisted 
cotunneling through the CQD as well as for singlet-triplet 
transitions. 

The precise manner in which these statements are 
quantified will now be explained. After completing the 
RG procedure, one arrives at an effective (or renormal- 
ized) Hamiltonian H in a reduced energy scale D, 

H = Hdot + Hi ea d + H cotun , (16) 

where the effective dot Hamiltonian is reduced to 

Aeiv 

written in terms of diagonal Hubbard operators, 

X AA = | A X A I- (18) 

At this stage, the manifold {A} <G N contains only the 
renormalized low-energy states within the energy interval 
comparable with Tk (to be defined below). Some of these 
states may be quasi degenerate, with energy differences 
\E\ — -Ea'1 < T K - However, T K itself is a function of 
these energy distances (see, e.g.j§i2ii22), and all the levels, 
which influence Tk, should be retained in l|17|) . 



The effective cotunneling Hamiltonian acquires the 
form 

H cot = J2 ( J o" Q ' S • S " Q ' + E J " Q ' R " • S " Q ' ) • ( 19 ) 

aa' \ n J 

Here S is the spin operator of CQD in its ground state, 
the operators s QQ represent the spin states of lead elec- 
trons, 

kk' aa' 

where f is the vector of Pauli matrices. In the conven- 
tional Kondo effect the logarithmic divergent processes 
develop due to spin reversals given by the first term con- 
taining the operator S. In CQD possessing dynamical 
symmetry, all R- vectors are involved in Kondo tunneling. 
In the following sections we will show how these addi- 
tional processes are manifested in resonance Kondo tun- 
neling through CQD. Note that the elements of the ma- 
trix J are also subject to temperature dependent renor- 
malization J% a ' -> J% a '(D /T). 

The cotunneling Hamiltonian (|19|) is the natural gen- 
eralization of the conventional Kondo Hamiltonian Js • S 
for CQDs possessing dynamical symmetries. In many 
cases there are several dot spin 1 operators depending 
on which pair of electrons is "active" . In this pair, one 
electron sits in well c and the other one sits in some well 
a. The other N — 2 electrons are paired in singlet states. 
This scenario applies if N is even. The spin 1 operator 
for the active pair is denoted as S a . (In some sense, the 
need to specify which pair couples to S = 1 while all other 
pairs are coupled to S — is the analog of the senior- 
ity scheme in atomic and nuclear physics (see, e.g.^)). 
The cotunneling Hamiltonian for CQD contains exchange 
terms Jq° S a ■ s aa . Then, instead of a single exchange 
term (first term on the RHS of Eq. (|19|l 'l. one has a sum 
Ja a • s™ . Additional symmetry elements (finite 
rotations and reflections) turn the cotunneling Hamilto- 
nian even more complicated. In the following sections 
we will consider several examples of such CQDs. It is 
seen from l)19fl. that in the generic case, both spin and 
R-vectors may be the sources of anomalous Kondo res- 
onances. The contribution of these vectors depends on 
the hierarchy of the energy states in the manifold. In 
principle, it may happen that the main contribution to 
the Kondo tunneling is given not by the spin of the dot, 
but by one of the R-vectors. 

Thus, we arrive at the conclusion that the regular pro- 
cedure of reducing the full Hamiltonian of a quantum dot 
in junctions with metallic leads to an effective Hamilto- 
nian describing only spin degrees of freedom of this sys- 
tem reveals a rich dynamical symmetry of CQD. Strictly 
speaking, only an isolated QD with N = 1 is fully de- 
scribed by its spin 1/2 operator obeying SU{2) symme- 
try without dynamical degrees of freedom. Yet even the 
doubly occupied dot with N = 2 possesses the dynamical 
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symmetry of a spin rotator because its spin spectrum con- 
sists of a singlet ground state (S) and a triplet excitation 
(T). Therefore, an R- vector describing S/T transitions 
may be introduced, and the Kondo tunneling through a 
dot of this kind may involve spin excitation under definite 
physical conditions, e.g., in an external magnetic field 6 . 
A two-electron quantum dot under Coulomb blockade 
constitutes apparently the simplest non-trivial example 
of a nano-object with dynamical symmetry of a spin ro- 
tator possessing an 15*0(4) symmetry. 

Dynamical symmetries SO(n) of CQDs are described 
by non-compact semi-simple algebras24. This non- 
compactness implies that the corresponding algebra o n 
may be presented as a direct sum of subalgebras, e.g., 
04 = 03 ® 03. Therefore, the dynamical symmetry group 
may be represented as a direct product of two groups 
of lower rank. In case of spin rotator the product is 
50(4) = SU(2) <g> SU{2). Generators of these subgroups 
may be constructed from those of the original group. The 
50(4) group possesses a single R-operator R, and the di- 
rect product is realized by means of the transformation 



S + R S-R 

2 ' 2 V ' 

Both vectors K and N generate SU(2) symmetry and 
may be treated as fictitious S=l/2 spins 21 . In some sit- 
uations these vectors are real spins localized in different 
valleys of CQD. In particular, the transformation (I21|) 
maps a single site Kondo problem for a DQD possessing 
50(4) symmetry to a two-site Kondo problem for spin 
1/2 centers with an SU{2) symmetry (see discussion in 
Refs. Il2lll3ri . For groups of higher dimensionality (n > 4) 
one can use many different ways of factorization, which 
may be represented by means of different Young tableaux 
(see Appendix D). 

Even in the case n = 4, the transformation l|21|) is not 
the only possible "two-spin" representation. An alter- 
native representation is realized in an external magnetic 
fieldi^. When the ground state of S/T manifold is a 
singlet (the energy 6 = Et — E$ > 0), the Zeeman split- 
ting energy of a triplet in an external magnetic field may 
exactly compensate the exchange splitting 5. This oc- 
casional degeneracy is described by the pseudospin 1/2 
formed by the singlet and the up projection of spin 1 
triplet. Two other projections of the triplet form the 
second pseudospin 1/2. The Kondo effect induced by ex- 
ternal magnetic field observed in several nano-objects 2 ^, 
was the first experimental manifestation of dynamical 
symmetry in quantum dots. 

We outlined in this section the novel features which 
appear in effective Kondo Hamiltonians due to the dy- 
namical symmetry of CQD exhibiting Kondo tunneling. 
In the following sections we will see how the additional 
terms in the Hamiltonian i|19|) influence the properties of 
Kondo resonance in various structures of CQDs. 



III. TRIPLE QUANTUM DOT IN PARALLEL 
GEOMETRY 

So far we have briefly mentioned a simple structure of 
CQD, i.e., double quantum dot with occupation N = 2 
and employed it to describe some generic properties of 
CQD enumerated in the previous section. This kind of 
an artificial molecule is the analog of a hydrogen molecule 
in the Heitler-London limitiSiH, and its 50(4) symme- 
try reflects the spin properties of ortho/parahydrogen. 
A much richer artificial object is a triple quantum dot 
(TQD), which can be considered as an analog of a linear 
molecule RH2. The central (c) dot is assumed to have a 
smaller radius (and, hence, larger capacitive energy Q c ) 
than the left (I) and right (r) dots, i.e., Q c ^> Qi, r . Fig. 
illustrates this configuration in a parallel geometry, where 
the " left-right" (I — r) reflection plane of the TQD is per- 
pendicular to the "source-drain" (s — d) reflection plane 
of metallic electrodes. 




FIG. 1: Triple quantum dot in parallel geometry and energy 
levels of each dot e a — e a — v ga (bare energy minus gate 
voltage) . 

To regulate the occupation of TQD as a whole and 
its constituents in particular, there is a couple of gates 
v g i,v gr applied to the I, r dots. The energy levels of 
single- and two-electron states in each one of the three 
constituent dots are shown in the lower panel of Fig. 1. 
Here the gate voltages v g i_ r are applied in such a way that 
the one-electron level e c of a c-dot is essentially deeper 
than those of the Z,r-dots, so that the condition H13fl is 
satisfied for the c dot, whereas the inequalities (|14|> and 
<P~5f) are satisfied for the "active" I and r dots. Tunneling 
between the side dots I, r and the central one c with am- 
plitudes Wi t r determines the low energy spin spectrum of 
the isolated TQD once its occupation N is given. This 
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system enables the exposure of much richer possibilities 
for additional degeneracy relative to the DQD setup men- 
tioned above due to the presence of two channels (I, r). 

The full diagonalization procedure of the Hamiltonian 
Hdot for the TQD is presented in Appendix E] When 
the condition (|15|l is valid, the low-energy manifold for 
N = 4 is composed of two singlets \Si), \S r ), two triplets 
\T a ) = \fJ-a) (a = I, r, \x a = l a , Q , l a ) and a charge trans- 
fer singlet exciton \Ex) with an electron removed from 
the c-well to the "outer" wells. Within the first order in 
(3 a <C 1 the corresponding energies are, 

E Sa = e c + e a + 2e- a + Q a -2W a /3 a , 

E Ta = e c + e a + 2e a + Q n , (22) 

E Ex = 2ei + 2e r + Q l +Q r + 2Wi0i+2W r p r , 

where the charge transfer energies in Ea. l|15|l (for deter- 
mining P a ) are E ac = Q a + e a — e c ; the notation a = l,r 
and a = r, I is used ubiquitously hereafter. 

The completely symmetric configuration, ei — e r = 
e, Ql = Qr = Q: Wi = W r = W, should be considered 
separately. In this case the singlet states form even and 
odd combinations in close analogy with the molecular 
states T,^ in axisymmetric molecules. The odd state S- 
and two triplet states are degenerate: 

E s+ = e c + 3e + Q - 4W(3, 

E s - = E Ta = e c + 3e + Q, (23) 

E Ex =4e + 2Q + iW0. 

Consideration of these two examples provide us with an 
opportunity to investigate the dynamical symmetry of 
CQD. 



A. Derivation and solution of scaling equations 

We commence with the case of TQD with even occu- 
pation N — 4 briefly discussed in Ref. ^3- This config- 
uration is a direct generalization of an asymmetric spin 
rotator, i.e., the double quantum dot in a side-bound 
geometrjsi^. Compared with the asymmetric DQD, this 
composite dot possesses one more symmetry element, i.e., 
the / — r permutation, which, as will be seen below, en- 
riches the dynamical properties of CQD. 

Following a glance at the energy level scheme 1(22(1 . 
one is tempted to conclude outright that for finite W, the 
ground state of this TQD configuration is a singlet and 
consequently there is no room for the Kondo effect to take 
place. A more attentive study of the tunneling problem, 
however, shows that tunneling between the TQD and the 
leads opens the way for a rich Kondo physics accompa- 
nied by numerous dynamical symmetries. 

Indeed, inspecting the expressions for the energy lev- 
els, one notices that the singlet states E$ a are modified 
due to inter-well tunneling, whereas the triplet states Ex a 
are left intact. This difference is due to the admixture 



of the singlet states with the charge transfer singlet exci- 
ton (see Appendix 0). As was mentioned in the previous 
section, the Kondo cotunneling in the perturbative weak 
coupling regime at T, e > Tk is excellently described 
within RG formalismi2i2&. According to general prescrip- 
tions of this theory, the renormalizable parameters of the 
effective low-energy Hamiltonian in a one-loop approxi- 
mation are the energy levels E\ and the effective indirect 
exchange vertices J^A' ■ 

To apply the RG procedure to the Kondo tunneling 
through TQD, let us first specify the terms Hi ea( i and 
H tun in the Anderson Hamiltonian (JSj. The most inter- 
esting for us are situations where the accidental degener- 
acy of spin states is realized. So we consider geometries 
where the device as a whole possesses either complete or 
slightly violated I — r axial symmetry. Then the quan- 
tum number a in Hi ea d l|10J> contains the lead index (s, d) 
and the channel index (l,r). If the axial I — r symmetry 
is perfect, the two tunneling channels are independent. 
To describe weak violation of this symmetry we intro- 
duce weak interchannel hybridization in the leads. This 
hybridization is characterized by a constant ti r <C Do, 
which is small first due to the angular symmetry, and 
second due to significant spatial separation between the 
two channels. The wave vector k is assumed to remain 
a good quantum number. Then, having in mind that in 
our model ekas = £kad = £fca, the generalized Hamilto- 
nian Ijl 0(1 acquires the form 

Hiead — ^ ^ ^ {^kan a bka + tlrC abkrJ Cdbka^ ■ (24) 
ka b—s,d a—l,r 

The tunneling Hamiltonian 1(121) is written as 

H tun = E(^>M^* AA + H - c -)- ( 25 ) 

AA ka ab 

We assume below V as = V a d = V a (see Fig0. 

The iteration processes, which characterize the two- 
step RG procedure contributing to these parameters are 
illustrated in Fig|2| The intermediate states in these di- 
agrams are the high-energy states \q) near the ultravio- 
let cut-off energy D of the band continuum in the leads 
(dashed lines) and the states |A) G — 1 from adjacent 
charge sectors, which are admixed with the low-energy 
states | A) G A by the tunneling Hamiltonian H t 1(12(1 
(full lines). For the sake of simplicity we confine our- 
selves with N = 3 states in the charge sector. 

In the upper panel, the diagrams contributing to the 
renormalization of Hdot are shown. In comparison with 
the original theorjiSS, this procedure not only results in 
renormalization of the energy levels but also an addi- 
tional hybridization of the states |A a ) via channel mix- 
ing terms in the Hamiltonian 1(25(1 . Due to the condition 
((13(1 . the central dot c remains "passive" throughout the 
RG procedure. 

The mathematical realization of the diagrams dis- 
played in Fig. 2a is encoded in the scaling equations for 
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A \ Ai / A 
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FIG. 2: RG diagrams for the energy levels Ea (a) and the 

effective exchange vertices J^a' (&) ( see text for further ex- 
planations) . 

the energy levels E\, 

^E A /dD = ^^-. (26) 

Here E A \ = E\ — E\, Ta are the tunnel coupling con- 
stants which are different for different A, 

TT a = 7rpo(V a 2 + 2V£), T Sa = a 2 a T Ta . (27) 

Here a a = yl — 2/3^, and po is the density of electron 
states in the leads, which is supposed to be energy in- 
dependent. These scaling equations should be solved at 
some initial conditions 

E A (D ) = E ( °\ (28) 

where the index (0) marks the bare values of the model 
parameters entering the Hamiltonian Ha ©• 

Besides, the diagram of Fig. 2a generates a new vertex 
M^ A , where the states A, A' are either two singlets Sj, S r 
or two triplets Tj,T r . The Haldane's iteration procedure 
results in a scaling equation 

dMir/d\nD = -7 (29) 

with the initial condition Mi r (Do) = and the flow rate 
7 = PaViV r tir. After performing the Haldane's procedure 
we formally come to the scaled dot Hamiltonian 

H = EA a X A « A " + Mi r X A " Aa (30) 

A„ A„A a 



with the parameters E\ a and M; r depending on the run- 
ning variable D. 

Due to the above mentioned dependence of tunneling 
rates on the index A, namely the possibility of Tt > Ts 
and Fs_ > Ts + , the scaling trajectories Ea{D) may cross 
at some value of the monotonically decreasing energy pa- 
rameter D. The nature of level crossing is predetermined 
by the initial conditions l|28|) and the ratios between the 
tunneling rates Ta- As long as the inequality \Ea\\ -C D 
is effective and all levels are non-degenerate, the scaling 
equations l|26(l may be approximated as 

TrdE A /d\nD = T A . (31) 

The scaling trajectories are determined by the scaling 
invariants for equations (|26|l . 

El = Ea(D) - tt-^a ln(^D/r A ), (32) 

tuned to satisfy the initial conditions. With decreasing 
energy scale D these trajectories flatten and become D- 
independent in the so called Schrieffer- Wolff (SW) limit, 
which is reached when the activation energies A a be- 
come comparable with D. The corresponding effective 
bandwidth is denoted as D (we suppose, for the sake of 
simplicity, that A a < Q a , so that only the states |A) with 
N' = N — 1 are relevant) . The simultaneous evolution of 
interchannel hybridization parameter is described by the 
solution of scaling equation (|29|l , 

M lr (D)= 7 \n^-. (33) 

If this remarkable level crossing occurs at D > D, we 
arrive at the situation where adding an indirect exchange 
interaction between the TQD and the leads changes the 
magnetic state of the TQD from singlet to triplet. Those 
states Ea , which remain close enough to the new ground 
state are involved in the Kondo tunneling. As a result, 
the TQD acquires a rich dynamical symmetry structure 
instead of the trivial symmetry of spin singlet predeter- 
mined by the initial energy level scheme H22|) . Appear- 
ance of the logarithmic enhancement of the hybridiza- 
tion parameter Mi r l|33|l does not radically influence the 
general picture, provided the flow trajectories cross far 
from the SW line , due to a very small hybridization 
7 < Ta < D. However, we are interested just in cases 
when the accidental degeneracy occurs at the SW line. 
Various possibilities of this degeneracy are considered be- 
low. 

The flow diagrams leading to a non-trivial dynami- 
cal symmetry of TQD with N = 4 are presented in 
Figs l3l5l6l The horizontal axis on these diagrams is the 
current energy scale D for lead electrons, where the ver- 
tical axes represent the energy levels E A (D). The dashed 
line E = —D establishes the SW boundary for these lev- 
els. 

Before turning to highly degenerate situations, where 
the system possesses specific SO(n) symmetry, it is in- 
structive to consider the general case, where all flow tra- 
jectories Ea{D) are involved in Kondo tunneling in the 
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SW limit. This happens when the whole octet of spin 
singlets and triplets forming the manifold (1221) remains 
within the energy interval ~ Tjc in the SW limit. The 
level repulsion effect does not prevent formation of such 
multiplet, provided ti r is small enough and the inequality 

M lr (D) < T K (34) 

is valid. At this stage, the SW procedure for con- 
structing the effective spin Hamiltonian in the subspace 
Ms = {Ti, Si, T r , S r } should be applied. This procedure 
excludes the charged states generated by H t to second 
order in perturbation theory (see, e.g.^). 

The effective cotunneling Hamiltonian can be derived 
using Schrieffer-Wolf procedure^! (see Appendix Q ■ To 
simplify the SW transformation, one should first ratio- 
nalize the tunneling matrix V in the Hamiltonian (|25|l . 
This 4x4 matrix is diagonalized in the s — d,l — r 
space by means of the transformation to even/odd com- 
binations of lead electron /c-states and similar symmet- 
ric/antisymmetric combinations of l,r electrons in the 
dots. The form of this transformation for symmetric 
TQD can be found in Appendix [BJ Like in the case 
of conventional QD 11 , this transformation eliminates the 
odd combination of s — d electron wave functions from 
tunneling Hamiltonian. 

It should be emphasized that this transformation does 
not exclude the odd component from H tU n in case of 
TQD in a series geometry^. The same is valid for the 
Hamiltonians (|24[l . (|25[) with ti r = 0: in this case the 
rotation in s — d space conformally maps the Hamiltonian 
Ha @ f° r TQD in parallel geometry onto that for TQD 
in series. Both these cases will be considered in Section 

Unlike the case of DQD studied in Refs. Il2ll3t where 
the spin operators are the total spin S and a single R- 
operator, describing S/T transitions, the TQD is repre- 
sented by several spin operators corresponding to differ- 
ent Young tableaux (see Appendix D). To order 0(|V| 2 ), 
then, 

H cot = Y,EA a X A ^ + ]T M lr X A ^ 

Aa A a Aq 

+ X! ( £ kanabka + UrC^ bka C abka ) 

ka b—s,d a—l,r 

+ ]T Jjs a ■ Sa + j lr P 

(35) 

a—l,r a— Z,r 

+ ^f T R-a ' S Q + Jl r ^ a ' Saa ' 

a—l,r a—l,r 

Here we recall that £a = E Aa (D), M lr = Mi r (D), and 
the effective exchange constants are 

T T _ v a rST _ T T 

J a ' J a ~ a a<J a , 

*_%("_!_ + _!_). (36) 
2 \e F - ei eF-£r/ 



The vector operators S a , R a , R a and the permutation op- 
erator P manifest the dynamical symmetry of TQD in a 
subspace Rs- The permutation operator 

P=J2( xSaSd + J2 (37) 

a=l,r /i=l,0,l 

commutes with S; + S r and R; + R,.. 

The spherical components of these vectors are defined 
via Hubbard operators connecting different states of the 
octet, 

S+ = V2(X^+X^), S~ = (fl+)t, 

s* = x uu ~ a: 1 " 1 - 

R+ = V2(X^-X S ^), R~ = (i?+) f , 

R z = _ {x o aSa + x s a o a) ^ 

R+ = V2{a- a X^ s - - a a X s ^), R~ = (fl+)t, 

Rl = ~{a- a X a ° s * +a a X s " *). (38) 

In addition to the spin operator l|20ll for conduction elec- 
trons, new spin operators are required, 

Saa = 2 C afc<T^ CT0 '' Cafe ' '' ■ $9) 

kk' ctct' 

An extra symmetry element (1-r permutation) results in 
more complicated algebra which involves new R-operator 
R and the permutation operator P interchanging I and 
r components of TQD. 

One can derive from the generic Hamiltonian l|35|) more 
symmetric effective Hamiltonians describing partly de- 
generate configurations illustrated by the flow diagrams 
of Figs. POSED These are the cases when the level cross- 
ing occurs in a nearest vicinity of the SW line in the flow 
diagram. It is important to distinguish between the cases 
of generic and accidental symmetry. In the former case 
the device possesses intrinsic I — r and s — d symmetry, 
i.e., the left and right dots are identical, the correspond- 
ing tunnel parameters are equal, and left and right leads 
also mirror each other, namely, e k i = £fer = e k- In the 
latter case the gate voltages violates I — r symmetry, e.g., 
make e; ^ e r , Vi ^ V r , etc. The level degeneracy is 
achieved due to competition between the I — r interdot 
tunneling and the lead-dot tunneling without changing 
the symmetry of the Hamiltonian. 

The basic spin Hamiltonian (|35|l acquires a more com- 
pact form, when a TQD possesses generic or acciden- 
tal degeneracy. In these cases the operators ()38|l form 
close algebras, which predetermine the dynamical sym- 
metry of Kondo tunneling. We start the discussion of 
the pertinent SO(n) symmetries with the most degen- 
erate configuration (Fig. |3J), where the TQD possesses 
generic I — r axial symmetry, i.e., the left and right dots 
are completely equivalent. Then the energy spectrum of 
an isolated TQD is given by Eas. (|23|l . The four-electron 
wave functions are calculated in Appendix l{A"|l . Such 
TQD is a straightforward generalization of the so called 
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FIG. 3: Scaling trajectories for P x £0(4) x 50(4) symmetry 
in the SW regime. Inset: Zoomed in avoided level crossing 
pattern near the SW line 



T-shaped DQD introduced in Refs. li3Tl29l It is clear, 
that attachment of a third dot simply adds one more el- 
ement to the symmetry group 50(4), namely the I — r 
permutation P, which is parity sensitive. 

To reduce the Hamiltonian l|35|l into a more symmet- 
ric form, we rewrite the Hubbard operators in terms of 
new eigenstates E\, recalculated with account of generic 
degeneracy (12. '{t and I — r mixing Mi r . In assuming 
that the latter coupling parameter is the smallest one, 
it results in insignificant additional remormalization ~ 
^f\Mi r \ 2 / (e + Q — e c ) of the states E$+ and Ee x - Besides, 
it intermixes the triplet states and changes their nomen- 
clature from left/right to even/odd. The corresponding 
energy levels are 



E T± {D) = E Ta TMir. 



(40) 



The flow trajectories for two pairs of states (T+,T_) 
and (S+,S~) diverge slowly with decreasing D. If this 
divergence is negligible in the scale of Tk, then three 
nearly coincident trajectories Et±,Es_ cross the fourth 
trajectory E$+ at some point, since the inequality Ts + < 
Tt± — Ts- with Tt± = STtpoV 2 , T$ + — oTt± is valid 

(a = yl — 4/3 2 ). If this level crossing happens near the 
SW line, we arrive at a case of complete degeneracy of 
the renormalized spectrum, and the whole octet Ms is 
involved in the dynamical symmetry (Fig. yj. The fine 
structure of the flow diagram in the region of avoided 
level crossing is shown in the inset. 

Since the tunneling occurs in even and odd channels 
independently, the parity is conserved also in indirect SW 
exchange. As a result, the effective spin Hamiltonian i|35|) 



acquires the form 

H cot = E K X K ^ 



£ kv c l ) ka C Vka 



A„ 



ka ri=g,u 



1=g,u ri=g,u 

? T \^S .« 4- W 7 ST R (1) -4- 7 ST R (2 h • « 
J 2 / j "1717 1 /_^\ J 2r) "'VV ' J in I r > 



(41) 



n 



Here tk g — £fc — U r , e/- u — £fc + Mr an d the lead operators 
c i-ika {v = 9i u ) a- re defined in Appendix IE1 The operators 
S r; , are defined analogously to S a , R a in Eq.lJSSJ, and 

the vector operators S™, Rl« , Rl^ are defined as: 



A w S r - 



R 



(i) 

7) 7) 



R 



(2) 



(42) 



The spherical components of the operators R^ and R^ 
are given by 



R 



(i)+ 



R 



"7777 \ it, ?7?7 / ) 

?( 2 )- = (p( 2 )+<|t 

= -x s ^\ R^ z = -X ^. (43) 

The spin operators for the electrons in the leads are 
introduced by the obvious relations 

1 



— 9 c lkcrTaa' Cgk' a' , 



kk' era' 



n E E C uka^ aa ' C uk'a' , 



(44) 



kk' 

ST ST J 



S S« ~~ 2 2 1 2. t C aka T w' c uk'o' , &ug — ( s gu) , 
kk' aa' 

instead of (|20() . Now the operator algebra is given by 
the closed system of commutation relations which is a 
generalization of the 04 algebra, 

[^I7J : ^rj'k] — i&jkm&rjr}' $r}m: 

[Rr/j : R-rj'k] — ^jkm ^7777' ^rjmi 

[Rrjj ? $r]' k] i^jkm&riri' Rrfm- (^^) 

The operators S,, are orthogonal to R r) , and the Casimir 
operators in this case are K n = S.2+R2 = 3. This justifies 
the qualification of such TQD as a double spin rotator 
which is obtained from the spin rotator considered in 
RefsiiSii^ by a mirror reflection. The symmetry of such 
TQD is P x 50(4) x 50(4). 

Four additional vertices appear in the effective spin 
Hamiltonian (|41|l at the second stage of Haldane- 
Anderson scaling procedure^. As a result, the exchange 
part of the Hamiltonian l|41|) takes the form 



Plant — 



r]=g,u 



»7=9,« 



r ST Tt -s 



+ r r Vs .s 4- Vf 7 ST R (1) 4- J ST Ti (2) ) ■ s 

-TJ 2 /_^ J VV S W ' / A J 2r) ""qrj ^ J 2fj ■"'jjjj J 

+ E J 3^ S ' ( ' s 'l + E J 3v R n ■ s n ■ ( 46 ) 
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The coupling constants in the Hamiltonian (|46|l are 
subject to renormalization. Their values at D = D are 
taken as a boundary conditions 



SO(4) 



JUD) = Jj(D) = Jg(D) = j£(D) 



V 2 



(47) 



Jl(D) = J^(D) = 0, jf g T {D) = aJf g {D) [i = 1,2,3) 

for solving the scaling equations. These can be written 
in the following form: 



df 



d\nd 



oi) 2 , (# T ^ 

1W T 2 + 2 



(ilnd 
dhid 

d\nd 

djil 

dhid 

aST 

dhid 



[( j3 T t) ) 2 + 2(-l)"m /l . J 3 T 1) 



CO' + 



o 2 T ) 2 , (.m?^ 



dj§ T 



Sr, j • 2 ■ 2 
[2j^% T + 2(-l)"m, r j^ + j 2 T jf, T 

2 [ J? (ifjf + .jfjf ) + 2 jf^ (ii^ + iifl ) 
»? 

[2j 3 T ^ T + 2(-l)"m ;rJ f, T + jfilf] , (48) 



where ji V — poJuj (i = 1,2,3), d = p^D and mi r ■ 
PqMi t . 

Solution of Eqs. (|48|l yields the Kondo temperature 



1 



T K = D 1 - 



8m/ r 



2m/ r 



(V3 + l)(3jT g +jf 5 T ) 



(49) 



The limiting value of this relation for perfect Z — r sym- 
metry of the leads is 

lim Tko = D cxp ( - - 4 — j . (50) 

Here and below the coupling constants j%{D) in all equa- 
tions for Tk are taken at D = D. We see that avoided 
crossing effect in the case of slightly violated I — r sym- 
metry of TQD turns the Kondo temperature to be a 
function of the level splitting H4U|) . Similar situation has 
been noticed in previous studies of DQDiSii^ and pla- 
nar QD with even occupation 6 , where Tk turned out 
to be a monotonically decreasing function of S/T split- 
ting energy 5 = Es — Et with a maximum at 8 = 0. 
Now the Kondo temperature is a function of two param- 
eters, Tk{Mij.,S). Looking at Fig|3] (which corresponds 
to 5 = 0) we notice that for large enough M; r , when 
the inequality 1|34|1 is violated, Mi r ^> Tk, the symmetry 




50(3) 



FIG. 4: Variation of Tk with parameters <5 and Mw (see text 
for further details). 

of TQD is reduced to 50(4) symmetry of S/T manifold 
with the Kondo temperature 



T K \ = D exp 



a ST 
Jl 



(51) 



Additional S/T splitting induced by the gate voltage 
(e; 7^ e r ) results in further decrease of Tk as a func- 
tion of S. The asymptotic form of the function Tk(S) is 



T K 
Tkx 



Tki 
5 



(52) 



(cfi&2i). In the limit of 5 — > D the singlet state should 
be excluded from the manifold, and the symmetry of the 
TQD with spin one in this case is SO (3). The general 
shape of Tk(Mi t ,8) surface is presented in Fig0] Thus 
the Kondo effect for the TQD with mirror symmetry is 
characterized by the stable infinite fixed point character- 
istic for the underscreened spin one dot, similar to that 
for DQDiZii 3 .. 

Now we turn to asymmetric configurations where E\ c ^ 
E rc , IY r 7^ Tt, ■ In this case the system loses the I — r 
symmetry, and it is more convenient to return to the 
initial variables used in the generic Hamiltonian (|35|) . 

When the Haldane renormalization results in acciden- 
tal degeneracy of two singlets and one triplet, Es l ~ 
E T , » E Sr < E Tr (Fig. EJ, the TQD acquires an SO (5) 
symmetry of a manifold {T/, Si,S r }. In this case the SW 
Hamiltonian l|35l) transforms into 

H= E A X AA + M lr (X s,Sr +X SrS <) (53) 

A=Ti,S t ,S r 

+ £kaC aka C aka + Ur ^ ^ C lke C a,k(j 

ka a— Z,r ka a=l,r 

+ JiS/ • S; + J 2 R; • S; + J 3 (Rl ' S r ; + R 2 ■ Sj r ), 

where J\ — jj , J 2 = Jf T and J3 = a,- J;,-. The spherical 
components of the vector operators Ri and R 2 are given 
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FIG. 5: Scaling trajectories resulting in an SO(5) symmetry 
in the SW regime. Inset: Zoomed in avoided level crossing 
pattern near the SW line. 



by the following expressions, 

R+ = -V2X Srh , R± = V2X Srh , R lz = -X Sr ° l , 

i? 2 + = (i?r) f , Rz=(Rfy, R 2z =Rl- (54) 

The group generators of the 05 algebra are the i-vectors 
S/ , R; from l|38f) and the operators intermixing I- and 
r-states, namely the vector R = Ri + R2, 

R+ = V2(X liS " - X sA ), R- = ( J R+) t , 

R z = -(X°' s " +X S " '), (55) 

and a scalar A interchanging Z, r variables of degenerate 
singlets 



A = i(X 



S r Si 



x s lSr) 



(56) 



The commutation relations Q},© in this particular 
case acquire the form 

[Slji •S'ifc] — i^jkmSlrrn [Rlj ; Rlk] = iS-jkmSlrrn 
[Rj,Slk] = iCjkmRmi [Rj > Rk] = i^jkmSlmi 

[Rlj, Slk] = iejkmRlrn, [Rlj,Rk]—iSjkA, (57) 

[Rj , A] = iRij , [A, R tj ] = iRj , [A, S tj ] = 0. 

The operators R/ and R are orthogonal to S; in ac- 
cordance with (JSJ). Besides, R; • R = 3X SlSr , and the 
Casimir operator is K. = Sf + R 2 + R 2 + A 2 = 4. 

Like in the case of double SO(A) symmetry studied 
above, the second step of RG procedure generates ad- 
ditional vertices in the exchange part of the interaction 
Hamiltonian l|53|l . 

Hcot = Jl$l ■ S; + J 2 R; • S/ + J 3 (Ri • S r ; + R2 ' Sir) 

+ J4S1 ■ S r + J 5 R • S; + Je(Rll • S r i + R-21 ' Sj r ) 

+ J7S1 ■ (Sir + S r ;) + Jg(Rl ' Si r + R2 ' S r i) 

+ JqR-I ■ S r + J10R ' S r + JllR/ ' (Sir + Sri) 

+ Ji 2 (Ri; • S/ r + R 2i ■ S,.;), (58) 



where R lt = X s ' Sr B.i, R 2 ; = R 2 X SrS '. The scaling 
properties of the system are determined by a system of 
12 scaling equations with initial conditions 



Ji(D) = Jj 



MD) = jf T 



'1 1 • j 2\j^i — 1 , (59) 
MD) = OrJ lr , MD) = (t = 4- 12) 

(see Eq. ISHlfor definitions) specifically 

djl r -2 , -2 , -2 1 -2 1 -2 1 • /„• 1 • \ 

= _ bl +^2 +J5 +J7 +J11U6 +J12) 



+ 



dj2 
dhid 

dhid 

dhid 



dfe 
dhid 

dje 
dhid 

dj7 
din d 

djs 
dhid 

djg 
dhid 
djw 
dhid 
dj 11 
dhid 
dj 12 
dhid 



-[2(iij 2 + Ma) + jrije + hi) - m irjs], 

-[j'sO'l + 3i) + JVO'5 + ho) ~ m lr(j6 + jll)}, 

- [ji + 37 + it + Jxo + jii + 3ii (ie + ji2 ) 



ji+3l+ 3l+3i2 



■[2jij 5 + j 7 (j 3 + js) - mirh], 
'[je(ji +34) - mtrja], 

(j3 + js)(j5 + 3lo) + (j2 + jd)(j6 + jl2) 



37(jl +34.) + ill (j2 +39)] , 

-\ja(ji + 34) + h(ji, + 3io) ~ n%ir(jn + ju)], 
-[2(j 4 J9 + jVjn) + iVO'e + 312) - mirjw], 
-[2j4iio + h(ja + js) - mirjd], 

— [inCii +3i) +hi32 +jd)]> 

- [ji2 (ji +34,)- mirja] ■ (60) 



From equations H60|l. one deduces the Kondo tempera- 
ture, 



T K2 = D[l- 



1 



mir 



jl + 32 + ^(jl+j2) 2 +2jl 



V2>, 



mir .(6i) 



Similarly to the previous case, this equation transforms 
into the usual exponential form when the I — r symmetry 
of the leads is perfect, 



lim T K2 =De h + h + V(h + jz) 2 + 2j'i . ( 62 ) 

mi r —>0 

Upon increasing m; r , the symmetry reduces from SO(5) 
to 50(4). The same happens at small m; r but with in- 
creasing §i — E$t — ■ In the latter case the energy Es t 
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is quenched, and at Si 3> Tki Eq. (|61|l transforms into 
T K = S^M-Ui^l) +hW]-_ 1 } (cf-l 4 )- On the other 
hand, upon decreasing S r = Et t — Es t the symmetry 
P x 50(4) x 50(4) is restored at S r < T K0 . The Kondo 
effect disappears when Si changes sign (the ground state 
becomes singlet). 

The next asymmetric configuration is illustrated by the 
flow diagram of FigHjJ 

E 0.2 0.4 0.6 0.8 1 D 



0.15 0.2 0.25 0.3 0.35 0.4 




FIG. 6: Scaling trajectories for SO(7) symmetry in the SW 
regime. Inset: Zoomed in avoided level crossing pattern near 
the SW line 

In this case, the manifold {7}, Si, T r } is involved in the 
dynamical symmetry of TQD. The relevant symmetry 
group is 5*0(7). It is generated by six vectors and three 
scalars. These are spin operators S a (a = l,r) and R- 
operator R; (see Eq. I38|) plus three vector operators Ri 
and three scalar operators Ai involving l — r permutation. 
Here are the expressions for the spherical components of 
these vectors via Hubbard operators, 

R+ = V2(X^°' +X 0iI "), R\ = Jf l!l " -X 1 " 1 ', 
R+ = V2(X ll °- + X - 1 '), R\ = X 1 - 11 -A 1 ' 1 -, (63) 
R+ = ^(A 1 "^ - X s ' 1 -), &l = -pf°- Si + X s '"-). 

The scalar operators A±, A^, A3 now involve the I — r 
permutations for the triplet states. They are defined as 

A 3 = z(X 0i0 "-A ^). (64) 

The (somewhat involved) commutation relations of oj 
algebra for these operators and various kinematic con- 
straints are presented in Appendix [D] The SW transfor- 



mation results in the effective cotunneling Hamiltonian 
H cot = £ E A X AA + Mi r (X T ^+X T ^) 

A=T(,S(,T r 

ka a—l,r ka a—l,r 

+ JlaSa ' S a + J2 S QQ ■ S aa 

a—l,r a—l,r 

+ J 3 (R^ 1} • s ri + Rf 5 • Sir) + J 4 R • Si, (65) 

where Ji„ = Jj, J 2 = Ji r , J3 = aiJ tr , J 4 = a t jf and 
S a a = J2fj. A PaAla S a - The spherical components of the 
vector operators Ri and R2 are 

Ri 1)+ = V~2X lrSl , fig*- = -V2X^ S \ 
R™+ = (R^y, i$>- = (R^ + )\ 
£« = -X°^, Rf] = (m^ (66) 

It is easy to see that S; r + S r ; = Ri + R2 and R3 = 
5(1) .5(2) 

« i + W 3 

Like in the case of 50(5) symmetry, the tunnel- 
ing terms Mi r X TaTii generate additional vertices in the 
renormalized Hamiltonian H cot . The number of these 
vertices and the corresponding scaling equations is too 
wide to be presented here. We leave the description of 
RG procedure for SO (7) group for the next section (as 
well as the case of TQD with odd occupation), where 
the case of M; r = is considered. In that case the scal- 
ing equations describing the Kondo physics of TQD with 
SO(n) symmetry are more compact. 

B. Section summary 

The basic physics for all SO(n) symmetries is the same, 
and we summarize it here. We have analyzed several ex- 
amples of TQD with even occupation in the parallel ge- 
ometry (Fig. [3 ■ Our analysis demonstrates the principal 
features of Kondo effect in CQD in comparison with the 
conventional SQD composed of a single well. These ex- 
amples teach us that in Kondo tunneling through CQD, 
not only the spin rotation but also the "Runge-Lenz" 
type operators R and R are involved. Physically, the 
operators R describe left-right transitions, and different 
Young schemes give different spin operators in the effec- 
tive co-tunneling Hamiltonians (see Appendix 

IV. TRIPLE QUANTUM DOT IN SERIES 

A. Motivation 

It was mentioned already in Section II I II that a TQD 
with leads possessing perfect I — r symmetry can be 
mapped onto a TQD in a series by means of geometrical 
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conformal transformation. Indeed, if the inter-channel 
tunneling amplitude t\ r in the Hamiltonian (|24[l is zero, 
one may apply the rotation in the source-drain space sep- 
arately to each channel and exclude the odd s — d com- 
bination of lead states both in the I- and r-channeli-i 
Since now each lead is coupled to its own reservoir, and 
one arrives at the series configuration shown in Fig [7] 

It is virtually impossible to conceive an additional 
transformation after which the odd combination of lead 
states are excluded from the tunneling Hamiltonian 28 . 
As a result, the challenging situation arises in case of 
odd occupation N = 3, where the net spin of TQD is 
S = 1/2, and the two leads play part of two channels 
in Kondo tunneling Hamiltonian. Unfortunately, despite 
the occurrence of two electron channels in the spin Hamil- 
tonian, the complete mapping on the two-channel Kondo 
problem is not attained because there is an additional 
cotunncling term J; r S • S; r + H.c. (sf r is determined by 
(|39[l ) which turns out to be relevant, and the two-channel 
fixed point cannot be reached. And yet, from the point of 
view of dynamical symmetry the series geometry offers a 
new perspective which we analyze in the present section 
for cases of even and odd occupation. 

B. Even occupation 

Consider then a TQD in series (FigEJ) with four elec- 
tron occupation N = 4. The Hamiltonian of the system 
can be written in the form, 

H = ^E Aa X A * A « + ^,E X X XX +J2 E ^ct ka c bka 

A a A ka b—s,d 

+ EEt(^ + V "Mk J* AA + H.c.]. (67) 

AA ka 

Here |A), |A) are the four- and three-electron eigenfunc- 
tions l|A5() and (|A9|I . respectively; E\, E\ are the four- 
and three-electron energy levels, respectively; X XA = 
|A)(A| are number changing dot Hubbard operators. The 
tunneling amplitudes V* A = V a (\\d a a\A) (a = h r ) de- 
pend explicitly on the respective 3 — 4 particle quantum 
numbers A, A. Note that direct tunneling through the 
TQD is suppressed due to electron level mismatch and 
Coulomb blockade, so that only cotunncling mechanism 
contributes to the current. 

The generic Hamiltonian i|35|) simplifies in this case to 

H cot = J2EA a X A ^+J2 E ^ b ct ka c bka 

A a ka b=l,r 

^ J a S Q ' s a 4~ Jlr-P ^ S a • Saa 
a—l^r a— Z,r 

+ E J a T ~&a ■ Sa + Jlr E ^ ' S ««> ( 68 ) 
a— Z,r a—l,r 

(the notation l,r is used for the electron states both in 
the leads and in the TQD). The antiferromagnetic cou- 
pling constants are defined by (|36|1 . The vectors S Q , R a 




FIG. 7: Triple quantum dot in series. Left (I) and right (r) 
dots are coupled by tunneling Wi, r to the central (c) dot and 
by tunneling Vi :T to the source (s) (left) and drain (d) (right) 
leads. 

and R a are the dot operators l|38() , P is the permutation 
operator l|37l) . and the components of the vectors s a , s aa 
are determined in Eqs. (|20|) (with a = a = l,r) and 
l]39p. The vector operators S a , R a , R a and the permu- 
tation operator P manifest the dynamical symmetry of 
the TQD. 

We now discuss possible realization of P x 5*0(4) x 
SO(4), 50(5) and £0(7) symmetries arising in the TQD 
with N = 4. Due to the absence of interchannel mixing, 
the avoided crossing effect does not arise in the series 
geometry. As a result, all cases of high symmetry are 
characterized by the same flow diagram of Figs. 13151 and 
EJbut without avoided crossing effects shown in the insets. 

Let us commence the analysis of the Kondo effect in 
the series geometry with the case P x 50(4) x 50(4) 
where E\ c = E rc and Tx r — r*r, (FigEJ- In this case 
the exchange part of the Hamiltonian (|68|1 is a simpli- 
fied version of the Hamiltonian l|46|l with the boundary 
conditions l|47|) . The scaling equations are the same as 
(|48|l with mir = 0. Solving them one gets Eq.J5U|) for the 
Kondo temperature. 

When E Sl » E Tl « E Sr < E Tr (FigEJ, the TQD pos- 
sesses the SO(5) symmetry. In this case the interaction 
Hamiltonian has the form 

H cot = JiSi ■ si + J 2 R • s/ + J 3 (Ri • s r t + R 2 • si r ), (69) 

which is the same as in Eq. I|53|) with Ri , R 2 determined 
by Ea. (|54|l . Respectively, the effective Hamiltonian for 
the Anderson scaling is a reduced version of the Hamil- 
tonian (|58|l 

H co t = Jl^l ■ S; + J 2 Rj • Si 

+ J 3 (Rl • S W +R 2 • Sir) + J4S1 -B r , (70) 

with the boundary conditions l|59|) for Ji, i = 1 — 4. 
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The scaling equations have the form 
dji 



d\nd 

dj2 
dhid 

djz 
dhid 

dji 
dhid 



= -2jiJ 2 , 



-33(31 + .74), 



(71) 



Of course, Eas. (|71(l for the Kondo temperature yield the 
limiting value 1)6 2 |l. 

When E T , « £r r < E Sl ,E Sr (Fig®, the TQD pos- 
sesses the 50(7) symmetry. In this case the Anderson 
RG procedure adds three additional vertices in the ex- 
change part of the basic SW Hamiltonian l|65|) , 



H C ot — y Jla^a ' s a + J2 ^ ' S aa • S a g 
a—l.r a—l,r 

+ J 3 (Ri 1} • Brl + R { 3 } ■ S lr ) + J 4 R* • Sj 

+ J]] </5aS a ' S s + JqUi ■ S r . 



(72) 



The boundary conditions for solving the scaling equa- 
tions are 

J la (D) = jJ, J 2 (D) = Ji r , J 3 (D) = ai Ji r , (73) 
J 4 (5) = a x Jf, J 5a (D) = J 6 (D) = 0(a = l, r). 

The system of scaling equations 

dju 
d\nd 

dji r 



dhid 

dj2 
dhid 

dh 
dhid 

dlnd 
dju 
dhid 

djzr 
dhid 

dje 
dlnd 



3 2 U + 3 f+jl 
3i . 33 



j2ir T y T y 

J2 (jlf + jlr + jgj + 35r) + Mjj + jg) 

2 

[32(74 + Je) + 33(jlr + JSr)] , 



jii + J f+ji 



;2 
35r 



3l . j 3 1 



2 2 

[2J5/J6 + J2J3] 



(74) 



is now solvable analytically, and the Kondo temperature 
is, 



Tk = D exp 



2j+ + ^jl +3(j 2 +j 3 ) 2 j ' 



(75) 



where j + = + j 4 + j ir , j_ = j u + j A - j lr . 

Like in the cases considered above, the Kondo tem- 
perature and the dynamical symmetry itself depend on 
the level splitting. On quenching the 5; state (increas- 
ing Si r = Es l — Et t ), the pattern is changed into a 
P x SO (3) x 50(3) symmetry of two degenerate triplets 
with a mirror reflection axis. Changing the sign of 5i r 
one arrives at a singlet regime with Tk = 0. 
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FIG. 8: Phase diagram of TQD. The numerous dynamical 
symmetries of a TQD in the parallel geometry are presented 
in the plane of experimentally tunable parameters x = T; /Y r 
and y = Ei c /E rc . 

The results of calculations described in this subsec- 
tion are summarized in Fig[S] The central domain of 
size Tko describes the fully symmetric state where there 
is left-right symmetry. Other regimes of Kondo tunnel- 
ing correspond to lines or segments in the {x, y} plane. 
These lines correspond to cases of higher conductance 
(ZBA). On the other hand, at some regions, the TQD 
has a singlet ground state and the Kondo effect is ab- 
sent. These are marked by the vertically hatched do- 
main. Both the tunneling rates which enter the ratio x 
and the relative level positions which determine the pa- 
rameter y depend on the applied potentials, so the phase 
diagram presented in Fig[H] can be scanned experimen- 
tally by appropriate variations of V a and v ga . This is a 
rare occasion where an abstract concept like dynamical 
symmetry can be felt and tuned by experimentalists. The 
quantity that is measured in tunneling experiments is the 
zero-bias anomaly (ZBA) in tunnel conductance g& The 
ZBA peak is strongly temperature dependent, and this 
dependence is scaled by Tk- In particular, in a high tem- 
perature region T > Tk, where the scaling approach is 
valid, the conductance behaves as 



g(T) ~ m- 2 (T/TW). 



(76) 



As it has been demonstrated above, Tk in CQD is 
a non-universal quantity due to partial break-down of 
dynamical symmetry in these quantum dots. It has 
a maximum value in the point of highest symmetry 
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P x 5*0(4) x 50(4), and depends on the parameters 6 a 
in the less symmetric phases (see, e.g., Eqs. I5D1 1521 R)2l 
175(1 . Thus, scanning the phase diagram means chang- 
ing Tk{8o)- These changes are shown in Fig|5] which 



Tk/Tko 



xSO (4) 




SO (3) 



6 rl /T K 



FIG. 9: Variation of Kondo temperature with S r i = v gr — v g i. 
Increasing this parameter removes some of the degeneracy 
and either "breaks" or reduces the corresponding dynamical 
symmetry. 

illustrates the evolution of Tjc with S r i for x = 0.96, 
0.8 and 0.7 corresponding to a symmetry change from 
P x 50(4) x 50(4), 50(7) to P x 50(3) x 50(3) and 
from 50(5) to 50(4), respectively. It is clear that the 
conductance measured at given T should follow variation 
of Tk in accordance with l(76|l . 



Since Tq > T Bl ,T B2 , the scaling trajectories cross in a 
unique manner: This is the complete degenerate config- 
uration where all three phase trajectories E\ intersect 
[E Q (D*) = E Bl {D*) = E B2 {D*)} at the same point D*. 
This happens at bandwidth D = D* (Fig lTUjl whose value 
is estimated as 



where 



D* = D exp 



WfE rc + W^E 

Fjl c Hi rr . 



(80) 



W?E 2 +W 2 E? 



If this degenerate point occurs in the SW crossover re- 
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FIG. 10: Scaling trajectories resulting in 50(4) x SU(2) sym- 
metry of FQD with N = 3. 



C. Odd occupation 

We now turn our attention to investigation of the dy- 
namical symmetries of TQD in series with odd occupa- 
tion N = 3, whose low-energy spin multiplet contains 
two spin 1/2 doublets \B\^) and a spin quartet \Q), 

E Bl = £ c + £ t + e r - | [WiPi + W r /9 r ] , 

E B2 = e c + e, + s r - - [Wi(3i + W r (3 r ] , 

E Q = e c + ei+e r . (77) 

There are also four charge-transfer excitonic counterparts 
of the spin doublets separated by the charge transfer gaps 
~ Ei — e c + Qi and e r — e c + Q r from the above states 
(see Appendix^). 

Like in the four-electron case, the scaling equations 
(I31|l may be derived with different tunneling rates for 
different spin states (Tq for the quartet and T Bi {i = 1,2) 
for the doublets). 

Tq = irpo (V? + V r 2 ) , 

r Bl = 7l 2 r Q , r B2 = 72 2 r Q , (78) 

with 

7i = ^/l-|(A 2 +^f) 5 72 = \jl-\(fi + P 2 r )- (79) 



gion, i.e., if D* « D, the SW procedure involves all three 
spin states, and it results in the following cotunneling 
Hamiltonian 

H cot = Y,( J * S + J f T R)-s a! (81) 

a—L,r 

where S is the spin 1 operator and R is the R-operator de- 
scribing S/T transition similar to that for spin rotator—. 
The coupling constants are 

tT _ 4 7 1 |K | 2 ST „ jT / S9A 

3(e_F - £a) 

This is a somewhat unexpected situation where Kondo 
tunneling in a quantum dot with odd occupation demon- 
strates the exchange Hamiltonian of a quantum dot with 
even occupation. The reason for this scenario is the spe- 
cific structure of the wave function of TQD with N = 3. 
The corresponding wave functions |A) (see Appendix |A"|) 
are vector sums of states composed of a " passive" electron 
sitting in the central dot and singlet/triplet (S/T) two- 
electron states in the l,r dots. Constructing the eigen- 
states I A) using certain Young tableaux (see Appendix 
D), one concludes that the spin dynamics of such TQD 
is represented by the spin 1 operator S corresponding to 
the I — r triplet, the corresponding R-operator R and the 
spin 1 1 2 operator s c of a passive electron in the central 
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well. The latter does not enter the effective Hamilto- 
nian H cot l|81fl but influences the kinematic constraint 
via Casimir operator K. — S 2 + M 2 + s 2 = The dy- 
namical symmetry is therefore SO (A) x SU(2), and only 
the 5*0(4) subgroup is involved in Kondo tunneling. 
The scaling equations have the form, 



djla 

dlnd 

dj2a 

dlnd 



' [jla + J 2 a] 
-2jlaj2a, 



(83) 



where j la = p J* , j 2a = PoJa ( a = ^ r )- From E 1 s 
(1831) we obtain the Kondo temperature, 



T K = max{T K i,T Kr }, 



(84) 



with T Ka = D cxp [-l/(jia + ha)]- 

An additional dynamical symmetry arises in the case 
when D* > D. In this case the ground state of TQD 
is a quartet S=3/2, and we arrive at a standard under- 
screened Kondo effect for SU(2) quantum dot as an ul- 
timate limit of the above highly degenerate state. 



D. Section summary 

To conclude this section, it might be useful here to 
underscore the following points: 1) The difference be- 
tween series and parallel geometries of TQD coupled to 
the leads by two channels exists only at non-zero inter- 
channel mixing in the leads, U r ^ 0. 2) One may con- 
trol the dynamical symmetry of Kondo tunneling through 
TQD by varying the gate voltage and/or lead-dot tun- 
neling rate. 3) In the case of odd electron occupation 
[N = 3) when the ground-state of the isolated TQD is 
a doublet and higher spin excitations can be neglected, 
the effective low-energy Hamiltonian of a TQD in series 
manifests a two-channel Kondo problem albeit only in 
the weak coupling regime*^ To describe the flow diagram 
in this case, one should go beyond the one-loop approx- 
imation in RG flow equations^ 4) The nominal spin of 
CQD does not necessarily coincide with that involved 
in Kondo tunneling. A simple albeit striking realization 
of this scenario in this context is the case of TQD with 
N = 3, which manifests itself as a dot with integer or 
half-integer spin depending on gate voltages. 



V. ANISOTROPIC KONDO TUNNELING 
THROUGH TQD IN SERIES GEOMETRY 



more exotic situations where the effective spin Hamilto- 
nian is in fact a "Runge-Lenz" Hamiltonian in the sense 
that the vectors R alone are responsible for Kondo ef- 
fect. Actually, just this aspect of dynamical symmetry 
in Kondo tunneling was considered in the theoretical pa- 
pers cited in Ref. |y and observed experimentally in Ref. 
I2M in which the Kondo effect in planar and vertical QDs 
induced by external magnetic field B has been observed. 
In this section we lay down the theoretical basis for this 
somewhat unusual kind of Kondo effect. 

Consider again the case of TQD in series geometry with 
N = 4. In the previous sections the variation of spin sym- 
metry is due to the interplay of two contributions to indi- 
rect exchange coupling between the spins S Q . One source 
of such an exchange is tunneling within the CQD (ampli- 
tudes W a ) and another one is the tunneling between the 
dots and the leads (amplitudes V a ). An appropriate tun- 
ing of these two contributions results in occasional degen- 
eracy of spin states (elimination of exchange splitting), 
and various combinations of these occasional degenera- 
cies lead to the rich phase diagram presented in Fig. [S] 
A somewhat more crude approach, yet more compliant 
with experimental observation of such interplay is pro- 
vided by the Zeeman effect. This mechanism is effective 
for CQD which remains in a singlet ground state after all 
exchange renormalizations have taken place. The nega- 
tive exchange energy 8 a may then be compensated by the 
Zeeman splitting of the nearest triplet states, and Kondo 
effect arises once this compensation is completed From 
the point of view of dynamical symmetry, the degener- 
acy induced by magnetic field means realization of one 
possible subgroup of the non-compact group SO(n) (see 
Eq. [2] and corresponding discussion in Section [Hj . The 
transformation 50(4) — » SU{2) for DQD in magnetic 
field was discussed in Ref. fili . 



B. Quantum dot with SU(3) dynamical symmetry 

In similarity with DQD, the Kondo tunneling may be 
induced by external field B in the non-magnetic sector 
of the phase diagram of Fig. [H] A very peculiar Kondo 
tunneling is induced by an external magnetic field B in 
the non-magnetic sector of the phase diagram of Fig. |S| 
close to the 50(5) line. In this case, a remarkable sym- 
metry reduction occurs when the Zeeman splitting com- 
pensates negative 8i tT — Es ur ~ ■ Then we are left in 
the subspace of states {Tli, Si, S r }, and the interaction 
Hamiltonian has the form, 



A. Generalities 

In all examples of CQDs considered in the previous 
sections the co-tunneling problem is mapped on the spe- 
cific spin Hamiltonian where both S and R vectors are 
involved in resonance cotunneling. There are, however, 



H cot = (J 1 R z 1 + J 2 R z 2 )s z + ^J 3 i(R+s^+RYs+) 



(J 5 R Z + J 6 R z 2 )s z 



J4 (R 3 sf r + Ras z i) 
J 7 {Rt-s; + R[s+). (85) 
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Here 



Ji(-D) = J 2 (D) 



2jf 



J 3i (D) = jf T , 



J 3r (D) = a r J lr , Ji(D)=Q (i = 4-7). (86) 
The operators Ri, R2, R 3 and Ri are defined as, 



R{ = 




X SlSl ), R+ = X hSl , R[ = 




m = 


¥ xUU - 


X s ^), R+ = X 1 ^, i? 2 - = 




R 3 = 


2 


R 4 = ^-X S ^'. 


(87) 



We see that the anisotropic Kondo Hamiltonian (|85|l is 
quite unconventional. There are several different terms 
responsible for transverse and longitudinal exchange in- 
volving the R-operators which generate both S /T and 
S a /Sa transitions. 

The operators 1)8 7|) obey the following commutation 
relations, 



[Rlj,Rlk] = iejkmRlr, 

V3, 



[i?2j, i?2/c] — iGjkmR2m, 



[Rlj, i?2fe] — -^-(-^3 — Ri)5jk{l — Sjz) 



+ -ejkm\Rlm5kz + R2m$jz ^-<> mz (R 3 + Ri) 

1 \/3 



y/3 

[RljjRj] — -R4S3Z J-(R2x — iR2y)(3jx + i^jy), 

1 V3 

[R 2 j, R3] = ^R 3 S 3 z £-(Rlx + iRly){5jx + iSjy), 

1 V3 

[R2j,Ri\ = —-Rn5jz + —r-\Rix + iRly){$jx — i$jy), 



[R 3 , R4] — -(i?2 - R*)- 



These operators generate the algebra u 3 in the reduced 
spin space {Tl;, Si, S r } specified by the Casimir operator 

R-i + R| + R 3 + R4 = 2 ■ 

Therefore, in this case the TQD possesses SU(3) sym- 
metry. These R operators may be represented via the fa- 
miliar Gell-Mann matrices A^ (i — 1, 8) for the SU(3) 
group, 

R-t = ^(Ai + tAa), Ri = \ (Ai - i\ 2 ) , 



R z _ A 3 

Ri - y, 



RI = \(X 3 + V3X S ), 



1 

2 

V3 



R 2 = o( A 4 + «A 5 ), i? 2 =-(A 4 -iA 5 ), 



1 

2 

V3 



i? 3 = -^(A 6 + iA 7 ), Ri = — (A 6 - £A 7 ) ■ 



As far as the RG procedure for the "Runge-Lenz" ex- 
change Hamiltonian l|85Jl the poor-man scaling procedure 
is applicable also for the R operators. The scaling equa- 
tions have the form, 



dji 
dlnd 

dj 3 i 
dlnd 
dj 3r 
dlnd 

dji 
dlnd 

djh 
dlnd 

djr 
dlnd 



dh 
dlnd 



j 3 r 





V3. .1 


J3*(il + f) 


- — J3rj4 



hr{h + 2j 2 + J5 + 2j 6 ) - \/3j 4 (i3( + V277) 



/V3 
A 3 



\/2 
2 

$6 



J7 



J 3 n 



dlnd 

JTj3rj4 



(89) 



where j = We cannot demonstrate analytical so- 

lution of this system, but the numerical solution shows 
that stable infinite fixed point exists in this case like in 
all previous configurations. 

Another type of field induced Kondo effect is realized 
in the symmetric case of 6 = Es CJ — Et^ u < 0. Now the 
Zeeman splitting compensates negative 5. Then the two 
components of the triplets, namely Eti u cross with the 
singlet state energy Es g , and the symmetry group of the 
TQD in magnetic field is SU (3) as in the case considered 
above. 



C. Section summary 

It has been demonstrated that the loss of rotational in- 
variance in external magnetic field radically changes the 
dynamical symmetry of TQD. We considered here two ex- 
amples of symmetry reduction, namely 50(5) — > SU(3) 
and P x SO{A) x £0(4) -> SU(3). In all cases the Kondo 
exchange is anisotropic, which, of course, reflects the ax- 
ial anisotropy induced by the external field. These exam- 
ples as well as the SO (A) — > SU(2) reduction considered 
earlie r 12 ! 13 describe the magnetic field induced Kondo ef- 
fect owing to the dynamical symmetry of complex quan- 
tum dots. Similar reduction SO(n) — > SU(n') induced 
by magnetic field may arise also in more complicated 
configurations, and in particular in the parallel geome- 
try. The immense complexity of scaling procedure adds 
nothing new to the general pattern of the field induced 
anisotropy of Kondo tunneling, so we confine ourselves 
with these two examples . 

Although the anisotropic Kondo Hamiltonian was in- 
troduced formally at the early stage of Kondo physics 31 , 
it was rather difficult to perceive how such Hamiltonian 
is derivable from the generic Anderson-type Hamilto- 
nian. It was found that the effective anisotropy arises 
in cases where the pseudo-spin degrees of freedom (like 
a two-level system) are responsible for anomalous scat- 
tering. Another possibility is the introduction of mag- 
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netic anisotropy in the generic spin Hamiltonian due to 
spin-orbit interaction (see Ref. |32J for a review of such 
models). One should also mention the remarkable possi- 
bility of magnetic field induced anisotropic Kondo effect 
on a magnetic impurity in ferromagnetic rare-earth met- 
als with easy plane magnetic anisotropy^ This model 
is close to our model from the point of view of effec- 
tive spin Hamiltonian, but the sources of anisotropy are 
different in the two systems. In our case the interplay 
between singlet and triplet components of spin multi- 
plet is an eventual source both of Kondo effect itself and 
of its anisotropy in external magnetic field. Previously, 
the manifestation of SU (3) symmetry in anisotropic mag- 
netic systems were established in RefA It was shown, in 
particular, that this dynamical symmetry predetermines 
the properties of collective excitations in anisotropic 
Heisenberg ferromagnet. In the presence of single-ion 
anisotropy the relation between the Hubbard operators 
for S = 1 and Gell-Mann matrices A were established. It 
worth also mentioning in this context the SU (4) D SO (5) 
algebraic structure of superconducting and antiferromag- 
netic coherent states in cuprate High-T c materials^. 



VI. CONCLUSIONS 

We have analyzed the occurrence of dynamical sym- 
metries in complex quantum dots. These symmetries 
emerge when the dot is coupled with metallic electrodes 
under the conditions of strong Coulomb blockade and 
nearly degenerate low energy spin spectrum. It can be 
achieved either by an application of an external magnetic 
field or due to dot-lead tunneling which, as we have seen, 
results in level renormalization and the emergence of an 
additional symmetry. Although the main focus in this 
paper is related to the study of triple quantum dots, the 
generalization to other quantum dot structures is indeed 
straightforward. 

Since we were interested in a symmetry aspect of 
Kondo tunneling Hamiltonian, we restricted ourselves by 
derivation of RG flow equations and solving them for ob- 
taining the Kondo temperature. In all cases the TQDs 
possess strong coupling fixed point characteristic for spin 
1/2 and/or spin 1 case. We did not calculate the tunnel 
conductance in details, because it reproduces the main 
features of Kondo-type zero bias anomalies studied exten- 
sively by many authors (see, e.g.) 6 ' n i 12 i 13 i 14 ' 21 i 22 ). The 
novel feature is the possibility of changing Tk by scan- 
ning the phase diagram of Fig. [SJ Then the zero bias 
anomaly follows all symmetry crossovers induced by ex- 
perimentally tunable gate voltages and tunneling rates. 

The main message of our work is that symmetry en- 
ters the realm of mesoscopic physics in a rather non- 
trivial manner. Dynamical symmetry in this context is 
not just a geometrical concept but, rather, intimately re- 
lated with the physics of strong correlations and exchange 
interactions. The relation with other branches of physics 
makes it even more attractive. The groups SO(n) play 



an important role in Particle Physics as well as in model 
building for high temperature superconductivity (espe- 
cially SO{5)). The role of the group SU{3) in Particle 
Physics cannot be overestimated and its role in Nuclear 
Physics in relation with the interacting Boson model is 
well recognized. This paper extends the role of these Lie 
groups in Condensed Matter Physics. 
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APPENDIX A: DIAGONALIZATION OF THE 
DOT HAMILTONIAN 

Here we describe the diagonalization procedure of the 
Hamiltonian of the isolated TQDs occupied by four and 
three electrons. The dot Hamiltonian has the form, 

a— l,r,c (7 a 

+ J2 ( W A<rd aa + H.c). (Al) 

a— Z,r 

(a) Four electron occupation: 

The Hamiltonian (|A1|) can be diagonalized by using the 
basis of four-electron wave functions 

it a ,i) = 4 T 4 T 4 T 4jo), M) = 4AAAi\°)> 

K0) = ^(4r4i + 4i^)4Ai\°)> 

\sa) = ^(dlA a i-4Aar)4Ani\°h 

\ex) = 4AArAn\°^ ( A2 ) 

where a = l,r; I = r, r = I. The Coulomb interaction 
quenches the states with two electrons in the central dot 
and we do not take them into account. The states (|A2|) 
form a basis of two triplet and three singlet states. In this 
basis, the Hamiltonian IjAlfl is decomposed into triplet 
and singlet matrices, 

fl .=(U). («) 

and 

/ ii V2W t \ 
H s =\ e r V2W r , (A4) 
\V2Wt V2W r e ex J 

where e\ — e c + e; + 2e r + Q r , e r = e c + e r + 2e; + Qi, 
and e ex — 2e; + 2e r + Qi + Q r - We are interested in the 
limit /?„ < 1 (/3 a are defined by Eq.(15)). So the secular 
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matrix may be diagonalized in lowest order of perturba- 
tion theory in f3 a . The eigenfunctions corresponding to 
the energy levels (22) are 



|S ) = y/\ - 2(31 W) - V2l3 a \ex), 

\Ta) = |t„) , (A5) 

\Ex) = ^l-2/32-2/3 r 2 |ea;) + V2pi\ Sl ) + V2f3 r \s r ). 

In completely symmetric case, Si — e r = e, Qi — Q r = 
Q, Wi = W r = W, the eigenfunctions corresponding to 
the energies l{2*3*)) are 



\s+ 
\s. 



) = yi~^ |S ' ) + ' Sr) -2/3|ex), 



V2 



\si) - \s r 



(A6) 



\T±) 

\Ex) = y/l^40\ex} + V20(\st) + \s r )), 



V2 



where (3 = W/(s + Q - e c ). 
(b) Three electron occupation: 

In this case the Hamiltonian (|A1|) can be diagonalized by 
using the basis of three-electron wave functions 



V6 



are the following combinations, 

\Bi,o) = Ji\h,a) - —0i\b r ,a) + — /3 r |6 ; , a), 
s/2 v/2 

1-82,0-) = 72|&2,cr) - _/3;|6 rjC r) - _ / 3 r |& ijCr ) ) 

3 1 

\Q,s z ) = \q,s z ) , s z =±-,±-, 



(A9) 



16,(7) = 



|6 c ,a) = --^ (d+d+ - d+d+ J d+ |0), 

g,±0=d+ ± d+ ± d+|O), 

1 \ _ (d+ ± d+ ± d+ + dj ± dt T d+ ± + d+ F d+ ± d+ )|0) 

2/~ V3 
|6 ic , a) = d+d+d+ |0), |6 rc , a) = d+ d+ d+ |0), 

|6;,o) = d+rf+d+ |0), |6 r ,a) = d+d+d+ |0), (A7) 



where cr =T,|- The three-electron states |A) of the 
TQD are classified as a ground state doublet \Bi), low- 
lying doublet \B%) and quartet |Q) excitations, and four 
charge-transfer excitonic doublets -B ac and B a (a = l,r). 
In the framework of second order perturbation theory 
with respect to (3 a (15), the energy levels E\ axe 



\B ac ,cr) = \Jl- 0i\b ac ,a) -Pa\ba,<r), 
\B r ,a) = ^l-2pf-P^\b r ,a)+f3 r \b lc ,a) 

+^P l {V3\b 1 ,a) + \b 2 ,a)) , 
\B u a) = yl\-2ffi-fi\b l ,a)+p l \b r ,c, a) 

-^/3 r (V3|6 1 ,(7)-|6 2 ,a)), 
where 71 and 72 are determined by Ea. l|79[) . 



APPENDIX B: ROTATIONS IN THE 
SOURCE-DRAIN AND LEFT-RIGHT SPACE 

In the generic case, the transformation which diago- 
nalizcs the tunneling Hamiltonian l|25l) has the form 



(Bl) 



With U a = V as IV a , V a = V ad /V a ) V* = \V as \ 2 + \V ad \ 2 

(a = l,r). In a symmetric case V as — V a d = V, this 
transformation simplifies to 

C-aeka 2 ^ (Casfca" ~l~ Cadka) 7 

Caoka 2 ^ ( ^askcr ~\~ Cadka^) 7 (1^2) 

and only the even (e) combination survives in the tun- 
neling Hamiltonian 



/ Cleka \ 






-( 


C-reka 




\ broker ) 





in 


Vl 









( Clska \ 


Vl 


Ul 


















u r 


V r 




Crska 1 








—v r 


U r J 





(B3) 



So the odd combination (o) may be omitted. 

If, moreover, the whole system "TQD plus leads" pos- 
sesses I — r symmetry, £; = e r , the second rotation in 
I — r-space 



E 



Bi 



6( + e r 



Eb 2 — £ c + e; + £r 



+ W r /3 r ] , 




Eq = e c + ei + e r , 

E Bac = e c + 2e a + Q a ~ W a (3a, 

E Ba = e a + 2e a + Qa + W a p a + 2WaPa- (A8) 

The eigenfunctions corresponding to the energy levels 



1 



transforms Hi eac [ + H tun into 
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1 







/ Qekcr 


1 


1 





1) 


Creka 








1 




dla 








-1 




V d ra 



(B4) 



Hlead + H tun = ^ [efcr^rjfea + V {c\ k(J d V cr + H.C.) (B5) 
71 k a 

With e^g = €k — Ur, £ku = £k + Ur- 
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APPENDIX C: EFFECTIVE SPIN 
HAMILTONIAN 



lations of the 07 Lie algebra, 



The spin Hamiltonian of the TQD with N = 4 occupa- 
tion in series geometry (FigEJ is derived below. The sys- 
tem is described by the Hamiltonian I|(i7l) . The Schrieffer- 
Wolff transformation^ for the configuration of four elec- 
tron states of the TQD projects out three electron states 
I A) and maps the Hamiltonian (|67|l onto an effective spin 
Hamiltonian H acting in a subspace of four-electron con- 
figurations (A| . . . I A') , 



H = e lS He- iS = H 



where 

* = -££ 



(j\ m 



A A 



X Ax c aka + H.c. 



AA (k)a,a 



(C2) 



Here (k) stands for the electron or hole states whose en- 
ergies are secluded within a layer ±D around the Fermi 
level. £aa = E\(D) — E\(D) and the notation a = l,r 
is used. The effective Hamiltonian with three-electron 
states I A) frozen out can be obtained by retaining the 
terms to order 0(|F| 2 ) on the right-hand side of Ea. (|Cl(l . 
It has the following form, 

H = ^E A X AA + £ tkac + aka c aka 

A (k)a,a 

EEE J^' a X AA 'c+ ka c ak>a/ (C3) 

AA'A kk' aa' a—l.r 



AA'A kk'crcr' 



[Saj,S a 'k] — 26 jkm^ ~*aa' 'ami [R[j > Rlk] — ^jkmSlm: 

[Rij, Sik] = iejkmRim, [Rij,S rk ] — [i?3j, Sik] = 0, 

\R-3j i R3k\ — i^jkmSrm: [Rfij i Srk] — i&jkmRsrrn 

[Rij,Rik] = iejkmS rm (l — Sj Z ){l — S kz ) 
i 1 

+ -^C-jkrnSl m {5jz + S kz ) — -(Slj5 kz — Sl k Sj z ), 
[R 2 j, R2k] = iejkmSl m (l - 8 JZ ){1 - 8 kz ) 

i 1 

~^T~^6jkmSrm{$jz ~t~ &kz) 7)(Srj&kz S r k^jz), 



[Rlj, R 2 k] 



[R 3 j , R 



lk\ 



[Rlj, Rlk] 
[i?2j, Rlk] 



( 2^-jkrn(,S rm Sjz 4" S[ m 5kz*) 
+ 2 [Slj^kz ~ S r k5jz + (Slz — S rz )Sj z S kz ] , 

■ iejkmRlmi^- ~ Sj z 

S kz 

- — (1 - S jz )Rij, 

■ ie jkm Ri m (S jz + + ^(1 - 5j z )Rij, 



^jkmRzm (&kz 



6 



2 ) ifi 1 ~ 5 kz )R 3 k, 



ie jkm R 3m (l-6 kz -f) 



+ 5 -f(l-5 kz )Rz k , 



where 



tAA' 
J kk'a 



(V AA )*V XA ' 

\ aa J aa 



1 



1 



E\\ — e ka E\/\ — e k i a 



aa' _ (v£ A rv#' 



kk'lr 



1 



1 



E\\ — £fe; -Ea'a — £fe' 



(C4) 



The constraint J^a X aa = 1 is valid. Unlike the con- 
ventional case of doublet spin 1 /2 we have here an octet 
A = {A;,A r } = {Si,Ti,S r ,T r j, and the SW transfor- 
mation intermixes all these states. The effective spin 
Hamiltonian (JOS) to order 0(|V| 2 ) acquires the form of 
Eq.Jgg). 



i\/2 - - 
[Ai, Sij] — iA 2 Sj Z H — — (Ri x 5j X — RiySjy), 



-iAl8j z ^-(RlySjx + RlxSjy), 



i\/2 - ~ 

[Al, S r j] — iA 2 Sj z — (R 2x 8j x — R2ySjy), 

i\/2 ~ 

[A 2 , S r j] = -iA x b 2Z + —^-{R 2y S jx + R 2x 5j V ), 

[A3, Stj] = -iR 2j (l - S jz ), [A 3 , Srj] = iRij{l - S jz ), 



APPENDIX D: SO(7) SYMMETRY 



The operators S;, S r , R;, Ri, R2, R3 and Ai (i = 
1,2,3) (see Eas.f33). (l(33"l) . (|o"i1 % l obey the commutation re- 



22 



[Ai 


Rij] - 


— '—^-{RzxSjx - R3y8jy), 


[A 2 


Rij] = 


—^-(RaySjx + R 3x Sj y ), 


[M 


-Rail = 


—— — (Rlx^ix ~ RlyS-iv), 




R3j] 




[A 3 


Rij] = 


-iRaJjz, [A 3 ,R 3j ] = 


[M 


R\j] = 






Ri ;1 = 


2 \ i ~ J ly u jx 1 ^IxVjy}-, 


[A3 
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l[S r x& jx H - S r yS ~'jy) ) 


[Ai 
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c\ \ 1 jjj 1 y J y / ' 


[A 2 


Rij] = 


— — i^S r ySj x ~h S rx 8jy)i 


[A3 
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[Ai 


A 2 ] = 


-i(Si,+S„), [A U A 3 ] 


[S a j , Rfj,k\ - 


- r a ^ v R 4- n a ^ n A 



(Dl) 



Here j,k,m are Cartesian indices, a = l,r; fj,,v = 1,2; 
n = 1,2,3; r^, a" fe and 0^ km are the structural 



constants, r, 



Their non-zero components are: 
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The following relations hold, 

S a • R ; = S Q • R 3 = 0, A X A 3 = A 2 A 3 = 0, 



1S1 



R 2 = + 3X SlSl , R 2 = + 3X Sl - 

A\ + A\ + A% = X» m + X^ r . (D2) 

Therefore, the vector operators S;, S r , R;, Rj and scalar 
operators Ai (i — 1,2,3) generate the algebra 07 in a 
representation specified by the Casimir operator 

2 3 
S? + S?+R? + £B*-Rj + Ri + X^ a =6. (D3) 

i=l i=l 



APPENDIX E: YOUNG TABLEAUX 
CORRESPONDING TO VARIOUS SYMMETRIES 

A TQD with "passive" central dot and "active" side 
dots reminds an artificial atom with inner core and ex- 
ternal valence shell. The many-electron wave functions 
in this nanoobject may be symmetrized in various ways, 
so that each spin state of N electrons in the TQD is 
characterized by its own symmetrization scheme. One 
may illustrate these schemes by means of the conven- 
tional graphic presentation of the permutation symme- 
try of multi-electron system by Young tableau 36 . For 
instance, triplet state of two electrons which is symmet- 
ric with respect to the electron permutation is labeled 
by a row of two squares, whereas the singlet one which is 
antisymmetric with respect to the permutation is labeled 
by a column of two squares. Having this in mind we can 
represent the singlet and triplet four electron states of 
the TQD (|X5|) by the four tableaux shown in Fig|TT] The 
tableaux Si (S r ) and T/ (T r ) correspond to the singlet and 
triplet states in which the right (left) dot contains two 
electrons (grey column in Fig lll|) whereas electrons in 
the left (right) and central dots form singlet and triplet, 
respectively. 



EB 



a 1 



iV2, a 2 



xy 



a 2 



a ly 



iV2. 



FIG. 11: Young tableaux corresponding to the singlet (S a ) 
and triplet (T a ) four electron states of the TQD. The grey 
column denote two electrons in the same dot (right or left). 

The Young tableaux corresponding to various SO(n) 
symmetries discussed in Sec. Ill can be obtained by com- 
bining the appropriate tableaux (Fig. I12|) . The high- 
est possible symmetry P x SO (A) x 50(4) is represented 
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PxSO(4)xSO(4) 1 — I I — I 

Px SO(3) xSO(3) 1 — I I — I 

SO(7) 1 — I I — I 

SO(5) ' — I 

SO(4) 1 — I 

FIG. 12: Young tableaux corresponding to SO(n) symme- 
tries. 



by four tableaux 7], T r , Si and S r since all singlet and 
triplet states are degenerate in this case. The symmetry 
P x SO (3) x SO (3) occurs when two triplets 7] and T r are 
close in energy and these are represented by the couple of 
Young tableaux in the second line. Following this proce- 
dure, the SO (7) symmetry can be described in terms of 
two triplets 7j, T r diagrams and one singlet Si diagram. 
Moreover, SO (5) symmetry is represented by two singlet 
Si, S r diagrams and one triplet T; diagram and, finally, 
one triplet and one singlet tableaux correspond to the 
SO (4) symmetry. 
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